Dark Matter in Minimal Trinification 

E. Di Napoli^^\ D. George^^\ M. Hertzberg^^\ F. Metzler^^^ and E. Siegel^^^ 



(1) University of North Carolina, 

Department of Physics and Astronomy, 
CB 3255 Phillips Hall, 
Chapel Hill NC 27599-3255, USA 
edodinOphysics . unc . edu 

(3) MIT, Center for Theoretical Physics, 
77, Massachusetts Av., 
Cambridge MA 02139, USA 
mphertz@mit . edu 

(5) University of Wisconsin, 
Department of Physics, 
1150 University Avenue, 
Madison WI 53706, USA 
siegelOphys . uf 1 . edu 



(2) The School of Physics, 
University of Melbourne, 
Victoria, 3010 Australia 
d . georgeSphysics . unimelb . edu . au 

(4) Institut fiir Theoretische Teilchenphysik, 
Universitdt Karlsruhe, 
76128 Karlsruhe, Cermany 
metzlerSparticle . uni-karlsruhe . de 



ABSTRACT 

We study an example of Grand Unified Theory (GUT), known as trinifi- 
cation, which was first introduced in 1984 by S.Glashow. This model has the 
GUT gauge group as [S'?7(3)]'^ with a discrete Z3 to ensure the couplings are 
unified at the GUT scale. In this letter we consider this trinification model 
in its minimal formulation and investigate its robustness in the context of 
cosmology. In particular we show that for a large set of the parameter space 
the model doesn't seem to provide a Dark Matter candidate compatible with 
cosmological data. 



* Article based on research supported by the US Department of Energy (grant DE-FG02- 
06ER41418) 



1 Introduction 



In this letter we would like to present the preliminary and suggestive results 
of a more ambitious and extensive research project. We study an example of 
Grand Unified Theory (GUT) in the context of certain requirements dictated 
by cosmology. In other words wc require the model in examination to ad- 
dress questions like "Is there a Dark Matter (DM) candidate? How abundant 
is this at present?" or "Can we find successful mechanisms for Baryogene- 
sis and Reheating?" . These questions arise from the more general program 
of using present cosmological data to constrain the enormous proliferation 
of phcnomenological works describing physics beyond the Standard Model 
(SM). 

We are going to concentrate on one of such models, known as "trinifica- 
tion", introduced for the first time in 1984 by S. Glashow [1], successively 
studied in detail by Babu et al. [2] and more recently by Willenbrock et 
al. [3] [4]. In particular we consider the minimal formulation of the model 
(respect to how the SM is embedded in it), and focus on the question of 
Dark Matter. We show that the model does not have a stable DM candidate 
compatible with VIdm < 0.24. We find that if we adjust the parameters 
such that we have a stable candidate, there is far too much DM. On the 
other hand, making the candidate unstable confiicts with Big Bang Nucle- 
osynthesis (BBN) constraints. In section 2 we are going to briefly present the 
model, list its salient features and focus on its advantages and downsides. We 
present our results in section 3. We conclude in section 4 discussing possible 
improvements of the model that would circumvent its difficulty in providing 
a viable DM candidate. 

2 Trinification in a nutshell 

The name "trinification" comes from its gauge symmetry: a triple replica of 
SU (3) conventionally written as 5'[/c;(3) xS't//,(3) x SUr{3)x'Ij2. The discrete 
group Z3 guarantees the gauge couphngs of the single SU{3) factors are the 
same at the GUT scale. The SM embedding is obtained by identifying the 
SUc{3) with the QCD gauge symmetry while the electroweak gauge group 
emerges as a result of breaking the other two SU (3) factors. Each SM fermion 
generation is embedded in a 27 = (1, 3, 3) © (3, 1, 3) © (3, 3, 1) representation 
of the gauge group. In order to better understand the field content we can 



1 



re-express this representation in terms of SM quantum numbers 
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We immediately notice that each generation includes the usual lepton and 

quark doublets and singlets plus some additional fields: we get two additional 
lepton doublets E and E"^, two neutral singlets Ni and N2 and two quark 
singlets B and B'^. We will see that the doublets under the breaking of 
[SU{^)f' will acquire a heavy mass 0{Mgut)) at tree level and become 
Dirac fermions. The singlets, instead, remain light and require a more refined 
adjustment involving radiative corrections. 

In order to have unification of the SM couplings at the GUT scale we 
also need two copies of scalars both in the same $"(27) = $^(1,3,3) © 
$gc(3, 1, 3) © $g(3, 3, 1) a = 1, 2 representation the fermions are in. Let's 
concentrate only on the and write its field content in terms of SM quantum 
numbers 

ru = (1,2, ~) cp-e = (1,2, 1) = (1,2, 

5«, = (1,1,0) S^, = (1,1,1) Si, = (1,1,0) 

The breaking [SU(3)]^ SUc(S) x SUl(2) x Uy(1) is obtained giving vevs 
to some of the singlets Sf,. The most general choice being (S2, are electrically 
charged so they cannot assume a vev) 

< ^3^^ >= Vi < Slg >= V2 < Sl^ >= -ug, (1) 

with f's ~ 0(Mgut = 10^^ Gel/) [3]. In the same fashion the electroweak 
symmetry is broken giving vevs to the electrically neutral components of the 
doublets (f)f charged under the SUl('2). A very general choice is 

< (<j>i,)' >= m < (0i,)o >= 712 < (4e)' >= ns 

< (0},)° >= U, < (<f>l,)' >= U2, (2) 
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with Ui, rii ~ 0{Mew)- The other scalar fields and $y do not acquire 
any vev since they carry color charge and would break SUc{5)- They are 
generically very heavy due to radiative corrections to their mass terms and 
do not show up in the low energy spectrum of the theory. From now on we 
will assume their masses to be of the order of Mqut- 



2.1 A simple case 

In a simplified version of the the model, we set ni = n2 = ns = ^3 = 
but keep all other vevs t(2, fi, '^'2) non-zero. The qualitative results will 
be exactly equal to the more general case and suffice to illustrate our point. 
A linear combination of the and four of the iS^" are eaten by twelve of 
the gauge bosons"*^ that become heavy with masses proportional to the Vi. 
Fine-tuning the quartic couplings, it is possible to obtain at most 5 light 
(~ Mew) Higgs doublets. At the same time Yukawa terms give masses for 
the fermions at the tree level. In general, such terms can be built pairing 
two fermion doublets and a scalar singlet or one scalar doublet contracting 
indices with one fermion doublet plus one one fermion singlet. In the first 
case we end up with a heavy (mass ~ Mqut) Dirac fermion meanwhile in 
the second case the fermions are light (mass ~ Mew)- Limiting our analysis 
to one fermion generation we obtain at tree level 

^ \l gfvf + glvl vfiE ^ \j h\vl + h\vl 

mu = giU2 rud ~ giui ~ hiUi (3) 



rrii, — h\U2 mjvi = hiU2 



h(uiU2 

ruE 



where ha and ga are couplings associated with the Yukawa terms proportional 
to 4>1^. This spectrum has some positive qualitative features and a negative 
one. As we have already mentioned B, B'^ and E, E'^ pair up to become very 
heavy Dirac fermions. The up and down quarks as well as the electron get 
different light masses. Unfortunately A^2, Ni and v are also hght, with the 
additional inconvenience that the last two seem to pair up to form a light 
Dirac fermion. This is highly undesirable and can be fixed invoking radiative 
correction induced by cubic scalar couplings [3] [4]. Calculating the one loop 



^The gauge bosons are in the 24 = (1, 1, 8) © (1, 8, 1) © (8, 1, 1) adjoint representation 
of [5C/(3)]3 
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contribution to mass terms for light fermions and using a seesaw mechanism 
the spectrum for these hght neutral fields become 



2 7-1 "'1 "'2 / A\ 

^ 9qF — (4) 

where F < Mqut is a factor of pure one-loop origin. In the presence of other 
two fermion generations there are additional iVj's and they can mix up and 
appear in a sort of hierarchy where the lightest of all can be pushed as far 
down as 10^ GeV . 



2.2 Light and darkness of Trinification 

We end this section by briefly reviewing the positive and negative features 
of this GUT model. Thanks to the two sets of scalar fields (six weak dou- 
blets this model achieves unifications of the SM gauge couplings at a 
scale around 10^^ GeV . The mechanism is similar to that of generic SUSY 
SU(5) GUT (although this unifies at IQi^^GeV): the six Higgs contribution 
to the /^-function is equivalent to the contribution of the two SUSY Higgs 
doublets and their fermion superpartners. The only difference resides in the 
fact that in Trinification there aren't gauge boson superpartners that are 
ultimately responsible of moving the GUT scale upwards. 

Usually, in other non-SUSY unification models there are mass degenera- 
cies between quarks and leptons since they come in the same representation 
of the gauge group. In Trinification this is avoided since quarks and leptons 
are in different representations in which the 27 is decomposed. In conclusion 
SM masses for fermions can be arbitrarily adjusted through Yukawa cou- 
plings to fit the experimental values. Prom this point of view Trinification is 
not any more predictive than the SM itself. Masses for the scalars are not 
protected and receive one loop quadratic corrections requiring fine tuning for 
at least one Higgs light doublet in Trinification does not provide any 
mechanism to solve the so called hierarchy problem. 

Unlike SU(5), in trinification gauge bosons conserve Baryon number and 
so do not mediate proton decay: the proton can only decay through Yukawa 
interactions. An acceptable value for its lifetime is recovered without the 
need of fine-tuning the involved Yukawa couplings. 

Due to abundant number of scalar fields, it is possible that baryogenesis 
is achieved at GUT scales [6] or at electroweak scales through a first order 
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phase transition [7]. Heavy scalars may be important in some Inflationary 
scenario. Moreover light neutral singlets may play the role of sterile neutrinos 
and be used to invoke a Leptogenesis mechanism [5]. The punch line being 
that this model offers a wide variety of possible developments in the context 
of Cosmology. In the next section we will address one of these issues. We 
investigate the possibility that the lightest of the neutral singlets can function 
as a viable candidate for Dark Matter and try to give an indicative answer. 

3 A Dark Matter Candidate 

Let us indicate the lightest neutral singlet as A'^. There are two main re- 
quirements that needs to satisfy in order to be a possible Dark Matter 
candidate. First, its relative abundance ft^^ has to be less than or equal to 
0.24. Second, its decay time needs to be much longer than the age of the 
Universe. 

3.1 Relative abundance 

The relative abundance at present (assuming the particle is stable) is cal- 
culated using the Boltzmann equation. The form of its solution depends on 
the regime (relativistic or not) at which it is approximated. To check if A^^ 
was non-relativistic at freeze out we need to verify under which condition the 
inequality = ^ > 3 is valid. If we assume that is non-relativistic the 
value of Xf depends logarithmically on the annihilation cross section and the 
mass of Ny^ [8] . 

Xf c^\n(^0m8-^MpLmN^{av)^ (5) 

Since we are interested only in estimating this expression we consider the 
dominant contribution to the annihilation channel. This is given by the 
following tree level Feynman graph^ 

^The contribution coming from two singlets annihilating in a virtual heavy X boson 
is even more suppressed since it is proportional to ^2"° and ^ vfggauge- 
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(6) 



Assuming 
for the cross-section 



and 



are light Higgs we are lead to the s-wave expression 



{av) 



(7) 



Here h is a generic Yukawa coupling. Plugging this expression in ^ and 
assuming the reasonable range 10^ GeV < rn^^ < 10^^ GeV we arrive at the 
conclusion that ^ < 10^ GeV in order to have Xf >3. Since Vi ~ Mqut the 
inequality cannot be satisfied implying that A^^^ is highly relativistic when it 
freezes out. 

We then equate the annihilation rate of A^^ to the Hubble rate H to 
obtain the freeze-out temperature 



Mpi{av) 



10^^ ^h-' GeV. 



(8) 



This gives a present-day abundance that is rather insensitive to (av) and 
much too large, namely 



(9) 



3.2 Decay time 

The previous estimate is based on the assumption that A^^^ is stable over 
the course of the universe requiring that tnx ^ 10^° yr. However, this is 
not correct for most reasonable choices of parameters in this model. We find 
that it decays not too long after it freezes out, and tends to destroy Big Bang 
Nucleosynthesis (BBN) unless we adjust some parameters appropriately. In 
order to estimate the lifetime of A^^^, we will consider only the dominant 
contribution, as we did earlier for the cross section. The most favorable 
decay channel is given by the following 
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— — y 




(10) 



In this graph is one of the scalars carrying color charge with mass of the 
order of Mqut- In calculating the decay rate it is a very good approximation 
to consider the decay products massless since they are highly relativistic. In 
such an approximation the decay rate is 



r^v. - TTT^^ (11) 



with Qq being the Yukawa coupling associated with the heavy colored scalar. 
For m^o^ ~ 10^^ GeV and rriN^ ~ 10^ GeV we get the following estimate for 
the decay time 

r^^ ~ lO^o^-^ GeV-^ ~ 10%"^ s. (12) 

For a typical value of Qq = 0.1 the decay time is about a 600 years making 
A^^ rather unstable and completely absent at present time in the universe. 
Moreover its decay products are so highly relativistic and get produced in 
such abundance that they destroy any product of BBN. We can reverse the 
reasoning and put an upper bound on the Yukawa coupling that stabilizes 

TN^ = 6.3 X IQ-^gq^ yr > 10^° yr ^ gg^2x 10"^ (13) 

This fine-tuning of is not so unreasonable, but it may backfire on the 
one-loop corrections to the mass of the neutral singlets compromising the 
efficiency of the radiative seesaw mechanism. 



4 Conclusion and discussion 

We have seen that the hghtest of the neutral singlets in the minimal trini- 
fication model is greatly overabundant at freeze out with the consequence 
that it is completely ruled out by cosmological data. This may be overcome 
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if the particle is not very stable but then it ruins BBN and creates another 
big problem. As it stands this GUT model does not seem to be a complete 
model of particle physics since it does not withstand one of the most needed 
cosmological requirements: the existence of Dark Matter. 

There is room for improvement though. On one side we may introduce 
a Z2 symmetry like the one introduced in SUSY models and make more 
stable without the need of adjusting the Yukawas opportunely. Such a dis- 
crete symmetry would forbid the Yukawa terms responsible for the dominant 
decay channel and stabilize the singlet. On the other hand we have to admit 
that our calculation of the cross-section for the relative abundance is a little 
"naive" since it doesn't take into account possible mixing of NiS with the 
other neutral Weyl fermions of the model. For example within just one gener- 
ations of leptons we have Eq E E and Eq G E^ plus the two neutral singlets, 
A'^i and N2-, and the neutrino v. The mass matrix of these five neutral fields 
is far from diagonal and can provide some mixing between them. The net 
result is that the lightest among them may have some non-zero coupling with 
the boson greatly enhancing the cross section amplitude. This mechanism 
may help in reducing the abundance to more accepted values. 

Finally we would like to thank the organizers of the LXXXVI "Les Houches" 
summer school, the staff and all our colleagues at the school that made our 
stay enjoyable and this work possible. 
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